Binding energy corrections in positronium decays 
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Abstract 

Positronium annihilation amplitudes that are computed by assuming a factoriza- 
tion approximation with on-shell intermediate leptons, do not exhibit good analytical 
behavior. We propose an ansatz which allows to include binding energy corrections and 
obtain the correct analytical and gauge invariance behavior of these QED amplitudes. 
As a consequence of these non-perturbative corrections, the parapositronium and or- 
thopositronium decay rates receive corrections of order and a^, respectively. These 
new corrections for orthopositronium are relevant in view of a precise comparison be- 
tween recent theoretical and experimental developments. Implications are pointed out 
for analogous decays of quarkonia . 
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Parapositronium (p-Ps) and orthopositronium (o-Ps) are bound states of e^e whose 
lifetimes for 27 and 87 channels, respectively, have been measured with high precision: 

r(p-Ps^77) = 7090.9(1.7) /is-^[l] (1) 

7.0514(14) 

r(o-Ps^777) = { 7.0482(16) (2) 

7.0398(29) . 

The corresponding theoretical predictions which include perturbative QED corrections to a 
non-relativistic treatment of the bound state wave function, have been computed also with 
high accuracy [^, ^: 

rp-Ps^77 = ^ 1- 5-— - + 2a2ln- + 1.75 30 - - — In^ - 

2 \ 4 vr a Vtt/ Itt a) 

= 7989.50(2) fxs~^ (3) 

r(o-Ps^777) = ^ 1 - 10.28661(1)--— In- + 5o - - — In^ - 

97r y 7r3« Kit / Zir a J 

= (7.0382 + 5o 0.39 x 10"^) /is~^ , (4) 

where m denotes de mass of the electron and a is the fine structure constant. As it can 
be observed, predictions for parapositronium are in very good agreement with experiment, 
while the experimental results for orthopositronium reported in Refs. [0, |^ largely disagree 
with the theoretical expectation. Recent theoretical efforts have focused on a more complete 
evaluation of the non-logarithmic O(a^) perturbative corrections for orthopositronium, with 
the result Bo = 44.52(26) [0] or, equivalently, 

r(o-Ps ^ 777) = 7.039934(10) fxs'^ . (5) 

This result renders even closer the theoretical prediction to the experimental measurement 
ofRef. i. 

Although these achievements of perturbative QED for bound states are impressive, little 
is known about the effects of non-perturbative corrections. In the present paper we imple- 
ment a mechanism, imposed by analyticity and gauge invariance, which allows to introduce 
binding energy (BE) effects in the positronium decay amplitudes. These new corrections 
affect the 77 decay rate of parapositronium at the order a^, while corrections of order 
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are induced in the orthopositronium decay rate and in the single photon spectrum of p- 
Dm— > e^e~7 and o-Ps^ 777 decays (p-Dm (o-Dm) denotes the J=0(1) bound state of the 
dimuonium system Moreover, the single photon energy spectrum in the last two 

decays previously mentioned is softened with respect to results obtained when BE effects are 
neglected. Our implementation of these non-perturbative corrections relies on well grounded 
basis of QED such as the requirement of a correct behavior of the decay amplitude in terms 
of the photon energies (analyticity) and on electromagnetic gauge invariance. We briefly dis- 
cuss at the end, the possible implications for analogous decays of quarkonium states. Details 
of our calculations will be given elsewhere 

Let us first try to explain how problems related to analyticity of the positronium decay 
amplitudes appears in the usual approaches. Current calculations of the positronium rates 
assume a factorization approximation of the dynamics contained in the decay amplitude. The 
amplitudes for positronium decays into a given final state X is approximated as the product 
of the amplitude for the bound state annihilation (described by $0? the e^e~ wavefunction 
at the origin), times the annihilation amplitude of free leptons into the final state X. Under 
this assumption, the positronium decay rates can be written as 

T{FsC'^'S) - X) = ■ {vreiaie^e- - X))„,.^,^o , (6) 

where a{e~^e~ — > X) denotes the total cross section for e+e^ —>■ X and Vrei is the relative 
velocity of the constituents in their center of mass frame. 

Since the leptons in the intermediate stage are taken on their mass-shell, the annihilation 
amplitude for e^e~ wy is automatically gauge-invariant. However, in this approximation 
the positronium decay amplitudes do not exhibit the expected analytical behavior. Indeed, 
for definiteness let us consider the 87 decay mode of orthopositronium where the photon 
energy is kinematically allowed to vanish. The behavior of the e'^e^ 87 amplitude in 
terms of the photon momenta is driven in the soft-photon limit by the intermediate electron 
propagators as follows: 

^ ^ A + m) 

where k is the four- momentum of the electron that emits a photon of four- momentum /j. 
Note that in the soft-photon limit Zj — > 0, the positronium decay amplitude seems to diverge 
as . Actually, this is only apparent since selection rules cancels these infrared divergencies 
in the static limit (k = (m, 0,0,0)), and the amplitude starts indeed at order l^. This is 
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in contradiction with the fact that in the soft-photon hmit the amphtude must vanish 
since o-Ps 777 involve only neutral external bosons. 

We can try to cure this bad analytical behavior by realizing that electrons in the interme- 
diate stage are always off their mass-shell due to BE effects. Actually, the e+e" bound states 
involve two mass scales: the mass M of the bound state and the mass m of the constituent 
electrons. These masses differ by terms of order and are related through the binding 
energy Ebi^d, 

Ebind = M-2m = --ma^ . (8) 

Thus, the relevant momentum scale for intermediate electrons is determined by M/2. In 
this case, the lepton propagator involved in the e^e~ 87 decay amplitude becomes (A;^ = 
MV4): 

i ^ Ah- A + m) , , 

fz- U-m -2li • - 72 ' ^ ^ 

where 7^ = — A/P/A ^ M^a^/IG, takes into account the bound state nature of the e^e~ 
pair in the rescattering process. The presence of 7^ in the denominator would provide to 
the decay amplitude a better analytical behavior. Unfortunately, the gauge invariance of 
the e~^e~ —>■ 777 amplitude is spoiled due to the presence of the two mass scales M and m. 
Thus, a more sophisticated procedure is required to restore analyticity without destroying 
gauge invariance. 

A natural way to incorporate BE effects is by considering a model with a loop of virtual 
leptons for the decays of positronium at lowest order. For definiteness we consider the 
77-7 decay of the positronium state B{'^'^'^^S): J = 0(1) being the spin of the p-Ps(p-Dm) 
state, n = 2 for p-Ps or p-Dm, and n = 3 for o-Ps or o-Dm decays (the e~^e~'y mode of 
paradimuonium can be reached from p-Dm^ 77* where 7* is a virtual photon). In this 
model we need to introduce a coupling FbT to describe the Be^e~ vertex, where F = 75(/)) 
for para(ortho)-positronium decay, 77^ is the polarization four-vector of the B{^S) state 
and Fb is a form factor that describes the structure of the vertex. The evaluation of the 
Bi^'^^^S) X decay amplitude follows standard rules (see Fig. 1). 

Applying the Feynman rules to Fig. 1 the decay amplitude is given by (contraction with 
photon polarization vectors e^(/i) • ■ ■ must be understood): 

where T'^'^ " is a properly symmetrized amplitude for annihilation of virtual lepton pairs into 
77.7. It is important to emphasize that this amplitude is gauge-invariant if Fb is constant. 
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To clearly illustrate how this loop model reproduces, in the limit of zero BE, the lowest 
order amplitudes of the on-shell approximation, let us consider the following ansatz for Fb'. 

FB = ^C% 2., ■ + 7') (11) 



where 7 contains the BE, $0 = ^(y^m^/Sn is the ground state wavefunction of e~^e at the 
origin, and C = — is a normalization constant. 



Using well known representations of Dirac-delta functions [|Tl| we can check that in the 
limit of zero BE (7^0, m — M/2 0) we obtain: 



(g-f)2-m2 (g + f)2_^2 2 ' go'-(Vq' + "^' + f )' 

Thus, upon (trivial) integration of Eq. (10) one gets: 

A<'^-(i?-n7) = ^Tr{r(/^-M)r^-(/^ + M)} , (13) 

where T^'^ denotes the reduced vertex evaluated at g = 0. Note that the factors (/P ± M) 
play the role of projector operators external to the action of photon vertices. 

Setting in the rest frame of positronium and using 2m = M (and the condition Pai]" = 
for orthopositronium) we arrive at the well known results of the factorization approxima- 
tion, namely: 

^-■■■(5 n7) = -C$o^Tr |l±3rr-j . (14) 

Let us return to the issues concerning gauge-invariance in the context of the present 
model when BE is not neglected. Gauge invariance requires that the amplitude in Eq. (10) 
(contracted with photon polarizations) vanishes when ea{li) ha for any external photon. 
The amplitude for p-Ps^ 77 is always gauge- invariant, no matter the specific form of Fb- 
This follows from the fact that Lorentz covariance for the axial-77 vertex implies that the 
amplitude should corresponds to an effective operator ~ F^^F^'^ , which is automatically 
gauge-invariant. This is not the case for o-Ps— > 777 decays. In this case the amplitude of 
Eq. (10) is gauge invariant provided Fb remains the same under the shifts g — » g + of the 
integration variable for all photons. Since this is achieved only if Fb is constant, it means 
that other contributions should be added to the orthopositronium decay amplitude in this 
loop model in order to compensate for this lack of gauge invariance. 

In view of these difficulties, we propose and ansatz for positronium decay amplitudes that 
fulfills gauge invariance and analyticity simultaneously. Our recipe contains three steps: 
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(a) evaluate the expression r^'^'" in Eq. (13) for constituents masses M/2 (this modified 
expression will be denoted by r^^^g); this will ensure gauge invariance, (6) introduce BE 
effects by multiplying the amplitude M.{B(^'^~^^S) wy) {i.e. the amplitude of Eq. (13) 
where we replace m by M/2 in the argument of the Trace operator) by a factor: Ai for p-Ps 
(p-Dm) states and 11^=1 -^i) for o-Ps (o-Dm), where we have defined Ai = P.li/{P.li + 7^) 
(note that Ai = A2 for p-Ps^ 27). This will warrant the correct analytical properties of 
the amplitude in the soft-photon limit and, (c) obtain the decay rate by integration over the 
physical phase-space determined by the mass M. 

Let us first note that the proposed ansatz is fulfilled automatically for the p-Ps— > 77 
decay. The reduced vertex in Eq. (13) is given by (remember g = and = M^): 

r = i-ie) ^ — ^ -2 — ^ — 

(f-/i) -m^ 



2 '-J 4 



In the second line of Eq. (15) we can replace m — > M/2 in the numerator because terms 
proportional to m in T^'^ cancel when performing the trace in Eq. (13). Therefore, gauge 
invariance is preserved independently of m and we have: 

A^'^^(p-Ps 77) = M'^p-Ps ^ 77) = ^Tr {75(/^ - M)T'Z/,{4^ + M)}Ai . (16) 

This amplitude satisfies analyticity, as required. 

Secondly, the reduced vertex for orthopositronium decay T'^'^^ can be worked in the 
following way. In order to accomplish gauge invariance, we are forced to replace m by M/2 
for the mass of the constituents and simultaneously add a third analytical factor Ai to 
each of the six amplitudes contributing to y'^'^^. This gives rise to the final gauge-invariant 
and analytical amplitude for orthopositronium decay: 

= {^(^ - ^^^mU^+ ^)} n (17) 



1=1 



where rja represents the polarization four-vector of orthopositronium. In the soft photon 
limit (/j 0), Ai^'^^ vanishes as required. 

Up to now we have considered the effects of nonzero BE corrections in the dynamics 
of positronium decays as expressed in the decay amplitude. It is clear that the physical 
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phase space for these decays is determined by the masses of external particles, in particular 
the initial available energy M. Thus, we will study the effects of these non-perturbative 
corrections in observables associated to positronium decays, which contain the BE effects in 
the dynamics (amplitude) and the kinematics (phase-space). 

27 decay of peirapositronium. 

The decay amplitude for p-Ps^ 77 obtained from Eq. (16) can be expressed as: 

M'\B{^S) ^ 77) = 2Ce^$o Ai ■ (18) 

The corresponding rate for this two-body decay where photons fly apart with energy M/2 
in the rest frame of p-Ps, is given by (2P./i = M^): 



1 (p-Ps 77) = -a m x 



1 



2 M2 
1 . 



I + 27VM2 



where we have used the definition of 7^ as given after Eq. (9) and we have neglected correc- 
tions of 0{a^). Thus, (non-perturbative) BE corrections for 27 decays of parapositronium 
are negligible small and are beyond present experimental precision. 

para-dimuonium (p-Dm) decay into e+e~7 

The decay mechanism for this process is similar to the previous one, but where one virtual 
photon converts into a e+e" pair. Therefore the corresponding amplitude is: 

AT'VDm ^ e-e-.) = 2Ce^^„^^pf^A,i^^I!}^, (20) 

where li and r = p + p' — P — h denote, respectively, the four-momenta of the photon and 
the e'^e~ pair. 

The single photon spectrum in this case is given by (a = Aml/M^, x = 2E^/M): 
cir(p-Dm ^ e+e-7) Wa^^l I ~, , I , , 

Observe that this spectrum falls as x^ when x ^ due to non-zero BE effects, instead 
of the usual behavior (proportional to x) expected when these effects are neglected. It is 
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interesting to note that when x <S Eq. (21) behaves as the corresponding photon 

energy spectrum in 7r° — > e"*"e~7 decay, for a point-hke pion vertex. Thus, BE corrections 
affect the shape of the spectrum or, conversely, actually probes the structure of the bound 
state. Indeed, Eq. (21) can be seen as a result that extrapolates the spectrum between 

the e'^e^'-f decay of a pscudoscalar point particle {n^ case) and the standard bound state 
calculations (paradimuonium decay without BE corrections). 

A closed (but long) analytic expression for the decay rate can be obtained from inte- 
gration of Eq. (21). A useful approximation that takes into account leading non- vanishing 
corrections of 0(q!^) is: 

fir 2 

r(p-Dm — > e'*"e~7) = 



Gtt 



aPy + 0(a^) 



(22) 



where the function Fq = (4/3)Vl - a(a-4) +2 ln[(l + Vl - a)/(l-Vl - a)] fixes the lowest 
order rate. 

37 decays of orthopositronium 

The squared unpolarized amplitude for o-Ps— > 777 obtained from Eq. (17), in terms of 
dimensionless photon energy variables — lE^^jM {x\ + ^2 + = 2) is given by: 

2 



E 1-^(37) roc 



X2XZ J \ xi^a / V x\X2 



\ — xxY ^ (\ — x^y ^ (\ — x^i 



Xi 



x- + ^ 



(23) 



As in the previous case, the single photon spectrum will be softened by BE corrections due 
to the last factor in the squared amplitude. 

The phase space integration of Eq. (23) can be performed in analytic form and expressed 
in terms of dilogarithmic functions. However, it is more illustrative to express the decay rate 
in terms of an expansion in powers of (or 7^). If we express the decay rate in terms of 
the mass m of constituents, we obtain: 



2(ti^ - 9) 

r(o-Ps 37) a^m- ^ ' 



1 ^ 2- 
1 — -a 



. (24) 
97r L 4 J ^ ' 

Thus, BE corrections affects the decay rate of orthopositronium at order which are indeed 

relevant when confronted to accuracy of present experiments. A comparison of Eqs. (4) and 

(24) indicates that the net effects of BE corrections is to resize the coefficient appearing in 

nonlogarithmic 0{a^) corrections of orthopositronium, namely: 

So ^ So - — fa 44.52(26) - 12.34 . (25) 
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i.e. a non-negligible 28% correction at order a^. 

Before concluding, let us address some comments on possible implications of BE correc- 
tions in analogous decays of quarkonia. Notice that BE in quarkonia is not simply related 
to the masses of the constituent quarks and to the coupling as in Eq. (8) due to confine- 
ment. If we assume, however, that BE corrections do affect quarkonium decays in a similar 
form as in positronium, we can address some apparent conflicts in some of their inclusive 
hadronic and radiative decays. First, the photon spectrum measured in J/ip ^ gg'y decays 
seems to be softer than predicted by perturbative QCD ||T2[ which indicate possible large 
non-perturbative effects. This is precisely the effect induced by BE corrections in the single 
photon spectrum of o-Ps— > 777. Second, different inclusive rates in quarkonium decays as 



discussed in |T3| will, in general, get decreased by BE corrections. This may increase the 
relatively low values of the strong coupling constants extracted from ratios of experimental 
quarkonium rates |TB|. Finally, these effects will also manifest in the "14 % rule" observed in 
the ratio BR{jp{2S) — > X)/BR{J/ip — > X) for single photon mediated decays [|14|, because 
BE are different for these two radial excitations of charmonium. 

In this paper we have computed binding energy corrections to positronium decays in 
a specific ansatz where these non-perturbative effects are introduced as a necessity to ac- 
count for analyticity and gauge invariance of the corresponding QED amplitudes. The BE 
corrections (of order a^) to orthopositronium decay are indeed relevant in view of recent 
efforts to achieve a precise comparison of theory and experiment. When extended to the 
quarkonium sector, these BE corrections may contribute to solve apparent discrepancies 
observed between experimental data and perturbative QCD calculations of some inclusive 
rates of quarkonia. 
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